The OPTX Project II: Hard X-ray Luminosity Functions of 
Active Galactic Nuclei for z < 5 

B. Yencho^ A. J. Barger^-^'^^ l. Trouille^, L. M. Winter^ 
ABSTRACT 

We use the largest, most uniform, and most spectroscopically complete to 
faint X-ray flux limits Chandra sample to date to construct hard (2 — 8 keV) 
rest-frame X-ray luminosity functions (HXLFs) of spectroscopically identified 
active galactic nuclei (AGNs) to z ~ 5. In addition, we use a new 2 — 8 keV 
local sample selected by the very hard (14 — 195 keV) SWIFT 9-month Burst 
Alert Telescope (BAT) survey to construct the local 2 — 8 keV HXLF. We do 
maximum likelihood fits of the combined distant plus local sample (as well as 
of the distant sample alone) over the redshift intervals < 2; < 1.2, < z < 
3, and < 2; < 5 using a variety of analytic forms, which we compare with 
the HXLFs. We recommend using our luminosity dependent density evolution 
(LDDE) model fits of the combined distant plus local sample over < 2; < 3 for 
all the spectroscopically identified sources and for the broad-line AGNs. 

Subject headings: galaxies: active — galaxies: evolution — surveys — X-rays: 
galaxies 



1. Introduction 



Understanding the evolution of the population of active galactic nuclei (AGNs) by means 
of a luminosity function is an area of continued interest, particularly for determining the ac- 
cretion history of supermass i ve black holes (SMBH s ) and the relationship to the galaxy pop 



ulation (jSalucci et al.l Il999l : lYu fc Tremaind l2002l : iMarconi et al.l l2004l : iBarger et al. 



2005 
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Merloni &: Heing|2008l ). While AGN surveys have been compiled in a variety of wavelengths, 
hard X-ray surveys are the least biased by absorption effects and are therefore very well suited 
to constructing models of the intrinsic space density of AGNs. In recent years, X-ray lumi- 



nosit y functions have been constructed from both hard (lUeda et al.l l2003l : iLa 



ranca et al. 



2005h a nd soft (iMivaii et al.l 120001 : iHasinger et 3.1.11200511 X-rav sample s, while ICowie et al. 
( l2003l ). iBarger et al.l (120051 ) . and more recently [Silverman et al.l (120081 ) have utilized both, 
taking advantage at redshifts z > 3 of the depth and small /^-corrections of soft X-ray data 
in measuring hard X-ray energies. 



Barger et al.l (120051 ) fitted their unbinned data over the redshift interval < 2; < 1.2 
and found that, in this region, the data were consistent with a pure luminosity evolution 
(PLE) model, which exhibits the same general shape at each redshift but shows evolution of 
the characteristic luminosity "knee", L^,. They noted that PLE does not continue to higher 
redshifts. Indeed, it has been shown that a more complex model with both luminosity and 



model, fits the higher redshift data 


Ueda et al. 


2003; 


Hasineer et al. 


2005 


2005; 


Aird et al. 


2008; 


Ebrero et al. 


20081; 


Silverman et al. 


2008). 



Also of i nteres t is the relationship o f the broad-line AGNs to the full population. 
Staff en et al. J2003h and lUeda et al.l (120031 ) observed that the fraction of broad-line sources 
(relative to the full population) increases with X-ray luminosity such that they dominate the 



AGN population at 



(iBarger et al 



2007 



ligh luminosities. Thi s has since been coii f irmed by a number of groups 



20051; 



La Franca et al 



Delia Ceca et al 



optical luminosity (IHao et al. 



20081: 



20051: 



12005 



Akvlas et al.ll2006l : iTozzi et al.ll2006l : iGilh et al. 



Hasinge ] 120081) and a similar correlation has been seen with 



Simpsonll2005l ). Having better HXLFs for both the full 



and broad-line data, consistent to high redshifts, will aid in understanding this effect. 

In this paper we construct full and broad-line 2 — 8 keV HXLFs using a data set drawn 
from the three high ly spectroscopically complete X-ray surveys of our OPTX project (see 
Trouille et al. I l2008h — the Chandra Deep Field North (CDF-N), the Chandra Large- Area Syn- 
optic X-ray Survey (CLASXS), and the Chandra Lockman Area North Survey (CLANS) — as 
well as from the ~ 1 Ms Chandra Deep Field South (CDF-S) and from an ASCA survey. 



While iBarger et all (l2005[ ) previously used the CDF-N, CDF-S, CLASXS, and ASCA data 
to construct their HXLFs, the OPTX project has since improved the spectroscopic complete- 
ness of the CDF-N and CLASXS fields and has introduced the CLANS field. The latter is a 
series of nine ~ 70 ks Chandra exposures covering a roughly ~ 0.6 deg ^ area consisting of a 
total of 761 sources, 533 of which have been spectroscopically observed (ITrouille et al.ll2008l ). 
With all the fields together w e have a sana p le of spectroscopically observed sources which 
is ~ 44% larger than that of IBarger et al. (l2005h . The set of spectroscopically identified 
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sources is similarly enlarged. With these additional data we can more accurately determine 
the binned luminosity functions and examine various model fits using maximum likelihood 
estimates. We also extrapolate our model fits to compare with a local 2 — 8 keV luminosity 
function that we construct from the first sensitive all-sky very hard (14 — 19 5 keV) X-ray 



survey in 28 years, the SWIFT 9- month Burst Alert Telescope (BAT) survey (ITueller et al. 



20081 : IWinter et al.ll2008l . l2009l ). We then combine the distant sample with the local sample 
and redo our model fits using the combined data set. The reason we do the model fits both 
with and without the local sample is because the BAT survey is not a 2 — 8 keV selected 
sample. However, at the 2 — 8 keV fluxes above ~ 2 x 10^^^ ergs cm"^ s^^ the BAT sur- 
vey should have chosen nearly all of the 2 — 8 keV sources over the whole sky, and we can 
therefore effectively construct a 2 — 8 keV sample. Comparison of the two sets of fits also 
gives insight into how reliable the extrapolations of the fits are outside the redshift range 
over which they have been made, since these complex multi-parameter fits may diverge from 
the observations quite rapidly outside the fitted region. 

In [|2] we briefly describe our distant X-ray sample and the local BAT sample that we 
use. In ^we construct rest-frame luminosities for both samples, and in §l]we determine the 
effective solid angles. We construct and analyze our binned luminosity functions in §3 and 
our unbinned maximum likelihood fits in ^ We summarize our results in ||71 

In this paper we use Hq = 70 km s~^ Mpc~^, Qm = 0.3, and = 0.7 and the convention 
that "log" is \ogiQ. 



2. X-ray Samples 
2.1. Distant Sample 



We use the 2 — 8 keV (hard) and 0.5 — 2 keV (soft) fluxes and redshifts of the X-ray 
samples from 5 highly spectroscopically complete surveys. These represent a wide range of 
area and depth regimes, with two very deep but narrow surveys (the CDF-N and CDF-S), one 
very large but shallow survey (ASCA), and two surveys of intermediate depth and coverage 
(CLASXS and CLANS). This complementary relationship allows for the construction of a 
HXLF which best describes the actual population of AGN. 

The PS 2 Ms CDF-N and ^ 1 Ms CDF-S (n ote that the latter field has recently been 
deepened by an additional 1 Ms; ILuo et al.ll2008l ) probe the deepest survey volumes over an 
area of ~ 0.1 deg^ each. Th e CDF-N flux limits ar e /2_8 keV ~ 1-4 x 10^^^ and fo.z 



1.5 X 10 ergs cm ^ s ^ (lAlexander et al.l 120031). and we take the CDF-S flux limits to 



1.5-2 keV 



be twice these values ( Alexander et al. 20031 : Giacconi et al. 2002). The optical spectra of 
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the X-ray sources i n the CDF-N were obtained by iBarger et a. 



Trouille et al.l (120081 ) . and those of the CDF- S were obtained by 



use the optical spectral classifications given in lTrouille et al.l (120081 ) and lSzokoly et al.l (120041 ) 



Szokolv et al.l (120041) . We 



(I200I I2OO3I . l2005h and 



to separate out the broad-line AGNs, which are defined therein as all sources which have 
optical lines with FWHM line widths in excess of 2000 km s~^. 



We include the ASCA sample of lAkiyama et al.l (120001 ) (hereafter referred to as simply 
"the ASCA sample") to provide us with higher fiux information in the 2 — 8 keV band 
over a large area of ~ 5 d eg^. The fiux limit of this 2—10 keV sample is /2- 



1.0 X 10 ergs cm ^ s ^ (lAkiyama et al.l I2OOOI ) . We convert the fiuxes and fiux limit of 
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this sample to the 2 — 8 keV band assuming a F = 1.8 spectrum. Of the 30 AGNs in the 
sample, only 5 are not classified as broad-line AGNs using our above definition. There is one 
ambiguous case in which the nominal FWHM of ~ 1800 km s~^ would fall below our limit, 
but this is based only on H/?, so we choose to keep this source as a broad-line AGN. 

We also include two intermediate-depth Chandra surveys of the Lockman Hole region. 
The first is CLA SXS, which is an ~ 0.4 deg^ field observed and analyzed by our team 
(lYang et al.ll2004l ). and the second is CLANS, which is an ~ 0.6 deg^ field initially observed 



by the Chandra /SW IRE team (P.I. B. Wilkes) and subsequently analyzed by our own team 
(iTrouille et al.l l2008l ). For both fields we adopt fiux limits of /2_8 keV ~ 3 x 10^^^ and 
/0.5-2 kcv ~ 5 X 10~^^ ergs cm~^ s"^ for our calcu lations. The optica l spe ctra of the X- 



ray so urces in the CLASXS field were obtained by ISteffen et al.l (1200 



(I2OO8I ) , and those of the CLANS field were obtained by 



the optical spectral classifications given in ITrouille et al 



) and Trouille et al. 



rouille et al.l (120081 ). We again use 



( 20081 ) to separate out the broad 



line AGNs. New and updated optical, near-inf rared, and mid - infrar ed data for the CLASXS, 
CLANS, and CDF-N surveys can be found in ITrouille et al.l (120081 ). 



In Table [T] we give a summary of the number and spectral classes of all the sources in 
each field. We note that the ASCA field was only observed in the hard band and that there 
are sources in the other fields which are measured in one band but fall below the detectable 
level of the other. Thus, in Table [2] we also give the number of sources of each spectral class 
in each field that have a valid fiux (measured and consistent with our chosen fiux limits) for 
the 2 — 8 keV and 0.5 — 2 keV energy bands separately. We hereafter refer to these as our 
hard and soft band samples. 

The completeness of each survey as a function of 2 — 8 keV and 0.5 — 2 keV fiux is given 
in Figures [T^ and [T]d, respectively. These values were determined by binning the sources 
according to fiux and finding the fraction of the spectroscopically observed sources which 
are spectroscopically identified. Due to the complicated observing program of the CDF-N, 
all of its X-ray sources must be considered as being spectroscopically observed. 



- 5 - 



While the completeness of each survey is quite good above 10~^^ ergs cm~^ s~^ in each 
energy band (70%-100%), there is a steady decline with decreasing flux, down to 50%-60% in 
some cases. This is a result of both the difficulty in obtaining spectra for faint sources and also 
the intrinsically lower fraction of broad-line AGNs at these fluxes, which are relatively easy 
to identify. Perhaps the most difficult sources to identify are the optically normal galaxies 
at z ~ 2 which occupy the intermediate-flux regime. This dip in completeness is apparent 
in the two deepest Chandra surveys. The corresponding increase in identiflcations at the 
faintest fluxes is associated with our ability to detect low-redshift star forming galaxies. 

To avoid problems associated with mixed identiflcation and classiflcation schemes, we 
choose to use only the spectroscopic identiflcations, being aware of our incompleteness at 
lower fluxes, rather than to supplement our spectroscopic i dentiflcations with p hotometric 



determinations or estimates based on hardness ratios (e.g., ISzokoly et al.l 120041 ) . This will 
manifest itself as an underestimate of the faint end of the luminosity function of the full 
population, particularly for high-redshift sources, but we will give estimates of the maximum 
possible error incurred from this approach. 
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Table 1. X-Ray Samples By Field And Spectral Class 



Category 


CDF-N 


CDF-S 


ASCA 


CLASXS 


CLANS 


Total 


503 


346 


32 


525 


761 


Observed 


460 


247 


32 


468 


533 


Identified 


312 


137 


31 


280 


336 


Broad-line 


39 


32 


25 


103 


126 


High-excitation 


42 


23 





57 


92 


Star-formers 


146 


55 


5 


77 


87 


Absorbers 


71 


20 





23 


22 


Stars 


14 


7 


1 


20 


9 



Table 2. Number of X-Ray Sources With Valid 2 — 8 keV/0.5 — 2 keV Flux Measurements 



Category 


CDF-N 


CDF-S 


ASCA 


CLASXS 


CLANS 


Observed 


299/413 


185/221 


31/0 


304/400 


378/419 


Identified 


184/288 


102/122 


30/0 


199/253 


267/283 


Broad-line 


37/38 


30/32 


24/0 


85/103 


109/123 


High-excitation 


33/39 


19/21 


0/0 


44/48 


78/74 


Star-formers 


74/134 


38/48 


5/0 


51/63 


62/59 


Absorbers 


36/63 


13/14 


0/0 


14/19 


15/18 


Stars 


4/14 


2/7 


1/0 


5/20 


3/9 
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2 - 8 keV Flux (ergs cm"^ s"') 




Fig. 1. — Fraction of observed sources that are spectroscopically identified for each survey in 
bins of (a) 2 — 8 keV flux and (b) 0.5 — 2 keV flux. Only bins which contain 10 or more sources 
are plotted (with the exception of the highest flux ASCA bin, which contains 8 sources but is 
included due to the survey's limited flux range). Due to the complicated observing program 
for the CDF-N, all sources are treated as spectroscopically observed. 
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2.2. Local Sample 



For our loc al comparison samp le we start with the 154 sources in th e SWIFT 9-month 
BAT sample of lTueller et~aD J2008h . for which Iwinter et all J2008I . boooh presented the X- 
ray properties, including the 2 — 10 keV fluxes, for 145. We then exclude the sources within 
the Galactic plane < 15°), where the optical identific ations are less complete. The 
remaining sky coverage is 30,500 deg^. iTueller et al.l (120081 ) compiled redshifts and optical 
spectral classifications for the sources. Here we classify all Seyfert Is through Seyfert 1.9s as 
broad-line AGNs. 

In Figure[2]we show the cumulative number counts per square degree versus the 2—8 keV 
flux (converted from the observed 2 — 10 keV flux assuming a photon index F = 1.8) for the 
|6| > 15° BAT sample [squares). We used the sky area covered by BAT at the 14 — 195 keV 
flux of each individual source from lTueller et al.l (120081 ). The number counts are consistent 
with the cumulative distribution of a uniform density of objects (—1.5 slope; red solid line) 
down to a 2 — 8 keV flux of 1.7 x 10^^^ ergs cm~^ s~^ [blue dashed line), which we take to be 
the completeness limit of the 2 — 8 keV sample selected by BAT. For our subsequent analysis 
we only consider the 2 — 8 keV sample above this limit, which we hereafter refer to as our 
local sample. All 37 of the sources in this sample have redshifts, 35 of which have z < 0.1. 
The median redshift for the sample is z = 0.016, and the mean redshift is z = 0.023. 26 of 
the 37 sources are broad-line AGNs. 
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10-'^ 10"" 10"'° 

F(2-8 keV) erg cm"^ s"' 

Fig. 2. — Cumulative number counts per square degree vs. 2 — 8 keV flux for the |6| > 
15° BAT sample of Tueller et al. (2008; squares) using the sky area covered by BAT at 
the 14 — 195 keV flux of each individual source. The red solid line shows the cumulative 
distribution of a uniform density of objects (i.e., —1.5 slope). The 2 — 8 keV sample selected 
by BAT is complete above about 1.7 x 10~^^ ergs cm~^ s~^ {blue dashed line). 
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3. Rest-Frame Luminosities 



3.1. Distant Sample 



To calculate the rest-frame 2 — 8 keV luminosities for our spectroscopically identified 
distant sample (see §2.11 for a description of the sample), a suitable redshift may be chosen 
such that the luminosities are calculated from the 0.5 — 2 keV (rather than the 2 — 8 keV 
band) fluxes for high redshifts. This is advantageous not only because the Chandra images 
are more sensitive at observed 0.5 — 2 keV, but also because at z = 3, observed 0.5 — 2 keV 
corresponds to rest-frame 2 — 8 keV. 

We call the redshift where we change over from using our hard band sample to using 
our soft band sample Zsput- Assuming a general photon index F = 1.8, we calculate the 
rest-frame 2 — 8 keV luminosities to be 



A J2/ \ I -^hard(l+2) if Z < 2;split , . . 

.. = «,(.)x^ ^^^,,[,,i^,)]-o. (1) 

where Fhard and Fsoft are the observed 2 — 8 keV and 0.5 — 2 keV fluxes in the hard and soft 
band samples, respectively. For our primary results we choose ^sput = 3.0. The K-corrections 
have different normalizations to characterize the redshift at which observed and rest-frame 
fluxes are equivalent. 

In Figure [3] we show rest-frame 2 — 8 keV luminosity versus redshift for the spectroscop- 
ically identifled sources in the distant sample, labeling each source according to the survey 
from which it was detected. The curves correspond to the flux limits for each survey, with 
the hard band sample flux limits being used for z < Zspm and the soft band sample flux limits 
being used for z > Zsput- This reflects our calculation of Lx and accounts for the resulting 
discontinuities. 

For comparison, in Figure H] we show the rest-frame 2 — 8 keV luminosities calculated 
from (a) the hard band and (b) the soft band samples separately, assuming a general photon 
index F = 1.8. The increased number of high- redshift sources in the soft band sample can 
be seen, though in both cases the samples are rather sparse at these redshifts. The soft 
band sample also has a greater number of low-redshift {z < 1), low-luminosity (Lx < 10^^ 
ergs s~^) sources relative to the hard band sample, but as we are in t erested in sources with 



Lx > 10^^ ergs s ^ (which we assume to be AGNs, IZezas et al.lll998l : iMoran et al.lll999l ). we 



will not be losing too much information by using only the hard band sample at low redshifts. 
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Redshift 

Fig. 3. — Rest-frame 2 — 8 keV luminosity vs. redshift for the spectroscopically identified 
sources in the X-ray surveys listed in the figure legend. The luminosities were calculated 
using the 2 — 8 keV fluxes from the hard band sample for z < 2;spiit and the 0.5 — 2 keV 
fluxes from the soft band sample for z > Zspia, with the i^'-corrections determined assuming 
a general photon index F = 1.8. The filled symbols designate broad-line AGNs, and the 
curves correspond to the flux limits of each survey (CDF-N - blue solid; CDF-S - purple 
long-dashed; both CLASXS and CLANS - red dot-dashed; ASCA - magenta short- dashed). 
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Fig. 4. — Rest-frame 2 — 8 keV luminosity vs. redshift for the spectroscopically identified 
sources in the X-ray surveys hsted in the figure legend. The luminosities were calculated 
using (a) the 2 — 8 keV fluxes from the hard band sample and (b) the 0.5 — 2 keV fluxes from 
the soft band sample, with the i^-corrections determined assuming a general photon index 
r = 1.8. The filled symbols designate broad-line AGNs, and the curves correspond to the 
flux limits of each survey (CDF-N - blue solid; CDF-S - purple long- dashed; both CLASXS 
and CLANS - red dot-dashed; ASCA - magenta short-dashed in panel a). 



- 13 - 



3.2. Local Sample 

For our local sample (see §2.21 for a description of the sample) we assumed a general 
photon index F = 1.8 and computed the rest-frame 2 — 8 keV luminosities from the observed 
2 — 8 keV fluxes, -Fhard, in this hard band sample using 

Lx = 47r4(^)Fhard(l + ^)~°-'. (2) 

In Figure [5] we show the rest-frame 2 — 8 keV luminosities versus redshift for our local 
sample {black squares) and for the spectroscopically identified sources in the ASCA {blue 
circles), CLASXS plus CLANS {red diamonds), and CDF-N plus CDF-S {green triangles) 
hard band samples. We show the luminosity limits corresponding to the flux limits of the 
surveys as the colored dashed lines. Our local sample is almost entirely a 2; < 0.1 sample 
(with two blazar exceptions, one of which is too luminous to appear on the plot) and is 
therefore essentially disjoint and complementary to our distant sample. The distant surveys 
are also complementary to one another and provide continuous coverage in both luminosity 
and redshift over the z < 5 redshift range where we are measuring the HXLFs. 
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0.001 0.010 0.100 1.000 10.000 

REDSHIFT 

Fig. 5. — Rest-frame 2—8 keV luminosity vs. redshift for our local sample {black squares; note 
that the blazar at z = 0.859 is too luminous to appear on the plot) and for the spectroscop- 
ically identified sources in the ASCA [blue circles), CLASXS plus CLANS {red diamonds) , 
and CDF-N plus CDF-S [green triangles) hard band samples. The luminosities were cal- 
culated using the 2 — 8 keV fluxes, with the i^'-corrections determined assuming a general 
photon index F = 1.8. The colored dashed lines show the luminosity limits corresponding 
to the flux limits of the surveys (from left to right: local sample; ASCA; CLASXS/CLANS; 
CDF-S; CDF-N). 
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4. Effective Solid Angles 

4.1. Distant Sample 

The sensitivity of X-ray detectors typically varies with the off-axis angle of the central 
pointing, resulting in an effective solid angle which is smaller for fainter sources. Our ap- 
proach to estimating the effective solid angle (as a function of flux) of our total hard and soft 
distant samples is as follows. We first bin our hard and soft samples separately according to 
(log) 2 — 8 keV and 0.5 — 2 keV flux, respectively. In each bin, we then compare the number 
of objects which are spectroscopically observed to independent determinations of expected 
number per unit area per unit fiux0. These differential number counts, here represented as 
n{F), use simple modeling of detector sensitivity to correct for the described bias. We there- 
fore estimate the effective solid angle for a bin with central value logF^ = 0.5(logFi-|-logF2) 
and observed count N as 



where X; = 2 — 8 keV or 0.5 — 2 keV and the differential number counts (in units of 
deg~^ per 10~^^ ergs cm^^ s^^) are given by 



nsiFs] 



{39 ±5){Fu/F, 
(32±2)(Fh/F, 

(12.49 ± 0.02) (Fs/Ffo 
(12.49 ±0. 02) (Fs/Ffo 



2.57±0.22 
1.63±0.05 



-2.5 

-1.7±0.02 



if Fh > Fb 
if Fh < Fh 

if Fs > F, 
if Fs < Ffe 



(4) 



for the = 2 - 8 keV fICowie et al.ll2002f ) and ^ = 0.5 - 2 keV bands (lYang et al.ll2004[ ). 



respectiveljl^. In each case, the break flux occurs at F), = 10 erg: 



s cm ^ s ^. 



^Note that for all of the fields except the CDF-N, the spectroscopic observations were targeted at the 
X-ray sources. Thus, we can assume that the unobserved sources would be like the observed sources and 
exclude them. In the method we are using to determine the area, when we remove the sources that were not 
observed, the area decreases accordingly. For the CDF-N, the observing strategy was more complex, and we 
need to treat all of the sources in that field as if they were observed. 



Yang et al.l (|2004l ) found it difficult to precisely fit the slope of their 0.5 — 2 keV data for fluxes > 



10 ergs cm ^ s "'^ but found it to be consistent with a fixed value of a = 2.5, resulting in our lack of 
errors on this value. 
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The resulting binned estimates are given in Figure [6], where the la Poissonian uncertain- 
ties are based on the number of objects in each bin. To reduce the dependence on binning, 
particularly for flux bins at the extremes of the distribution where the data are scarce, it is 
advantageous to fit the binned data with simple curves that show little variation with the 
details of the binning procedure. Smooth curves also facilitate the numerical computations 
necessary for estimating the luminosity function. Therefore we use the binned data to fit a 
truncated harmonic series in log Q — log F space according to a general linear least squares 
approach. We only assume these fits to be valid over the range of fluxes covered by the data. 
For fluxes above the brightest flux bin, we assume the effective solid angle is approximately 
constant and equal to the high flux endpoint of the fit; for fluxes lower than the faintest 
flux limit of our samples which are thus observationally inaccesible, we assign a value of zero 
to the curves. This process is performed separately for the hard and soft samples and the 
resulting curves are presented in Figured Calling these curves and r^g, we then define 
a single function which contains information from both distributions: 



n(\og Lx,z) 



Qsi^s) if z> Zspiit 



(5) 



Here Fh and Fs are the hard and soft fluxes that would be observed for a source with a 
rest-frame 2 — 8 keV luminosity Lx at redshift z, assuming a general photon index F = 1.8. 
This use of Zsput to define a single effective solid angle describing both the hard and soft 
samples follows from our definition of Lx- 

We also investigated how the use of different authors' fits to their differential number 
counts data would affect our res ults. Using the ri^LFw) an d ns(-Fs) de r ived f rom the 2 — 8 keV 
and 0.5 — 2 keV samples of both lTrouille et al.l (120081 ) and lKim et al.l (120071 ). we recalculated 
^h(-^h) and f2s(-^s) according to the above method. We find these forms of f2H(-FH) to be 
in very good agreement with our previous determination, but see some small discrepancies 
at both the faint and bright ends of fis(-Fs)- However, these differences had no significant 
effect on our final results, given the uncertainties. 



4.2. Local Sample 

Since the local sample is selected in the BAT high-energy band, determining the solid 
angle versus 2 — 8 keV flux relation is not straightforward. Above our chosen flux limit of 
1.7 X 10^^^ ergs cm~^ s~^ in the 2 — 8 keV band the situation is simpler, since most of the 
sources could be found across the full area of 30,500 deg^ for the sky at |6| > 15°. However, 
even at these 2 — 8 keV fluxes, about a third of the sources with z = — 0.1 (those with 
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higher flux ratios) have covered areas which are smaller than this, with the smallest being 
14,400 deg^. We have therefore tried two area versus 2 — 8 keV flux relations. In the first we 
assumed a constant area of 27,000 deg^ equal to the mean area in the sample. In the second, 
we derived an area versus 2 — 8 keV relation from the area versus 14 — 195 keV relation of 



Tueller et al.l (120081 ) by a simple conversion of the 14 — 195 keV flux to the corresponding 



2 — 8 keV flux found when assuming a fixed ratio of the two fluxes equal to the median ratio 
of 0.53. These fixed and flux dependent 2 — 8 keV area relations give almost identical results 
for the luminosity functions, and in the following we show the results for the fixed area case 
only. 
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Fig. 6. — The total effective solid angle of the distant surveys as a function of (a) 2 — 8 keV 
and (b) 0.5 — 2 keV flux. The points represent the binned data, with the error bars indicating 
the 1 cr Poissonian uncertainties. The solid curves are harmonic series fits to the data. 
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5. Hard X-ray Luminosity Functions: Binned 



5.1. IPage Carreral (120001) Method 



Let the hard X-ray luminosity function, defined as the number of objects per comoving 
volume per (log) luminosity, be given by 



0(log Lx,z) 



d^log Lx,z) 
dlogLx 



(6) 



We wish to estimate the value of this function in a bin of data defined by logLi < logL < 
logL2 and < z < z-?,. We begi n by following the derivation of the binned luminosity 
function by iPage fc Carreral (120001 ) . Given M different surveys, each with a unique flux 
limit and effective survey area, the expected number of objects in this bin would be given 
by 



M 



logL2 



logLi 



<{logL) 

0(log L, z)Qj (log L, z) — — (-z) dz dlog L 

az 



(7) 



where dV/ dz is the differential volume corresponding to a solid angle of 1 steradian, VLj (log L, z) 
^j{F) is the effective survey area of the j*^ survey, and z^^^(\ogL) is the maximum possible 
redshift within the bin for which an object of luminosity L could be detected by survey j 
(given the restrictions of that survey's flux limit). If we now assume that (p varies little over 
this bin, we can factor it out of the integral, yielding 



f-logL2 rz 
JlogLi Jz 



Qj (log L, z) {dV/dz) dz dlog L 



where is the observed number of objects in the bin. Wr itten in this form, the estimate is 
an implementation of the method of iPage fc Carreral (120001 ) which adheres to the "coherent" 
addition of samples described by lAvni fc Bahcalll (jl980l ). We can manipulate this further, 
though, to flnd a form more suitable for our data. 

Consider that each Qj (F) is a binned function derived according to the procedure given 
in ^ but using only the data of the j^^ survey. Then, by construction, each of these is equal 
to zero for any flux fainter than its flux limit: 



(9) 
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This means that the z{^^^{logL) provide redundant information, and we may set the upper 
hmit of the redshift integrals to Z2 without approximation. The sum can then be moved 
inside the integrah 



/S'C E,fi.(log^,^) (dy/dz)d^dlogL 



(10) 



The factor in brackets can be identified with the total effective survey area from Equation 
[5l The final form, then, for our estimate is 



N 



iLiL:i::^i^ogL,z)^iz)dzd\ogL 



dv, 



(11) 



with the 1 a Poissonian uncertainty (jGehrelslll986l ) on this value determined by the number of 
objects in the bin. This form is convenient because it packages all the information concerning 
the limitations of the surveys into a single quantity, Q{log L, z), rather than in the limits 
of multiple integrals. Also, given the binning required to determine the effective, flux- 
dependent solid angle for a particular survey, a direct calculation of the total solid angle for 
all surveys allows for the best possible statistics in each bin. This is the primary advantage 
of Equation [11] over Equation [HI 



5.2. Results for the Distant Sample 

We calculate binned HXLFs according to the above procedure for the distant sample 
having Lx > 10^^ ergs s~^ in three low-redshift bins and two high-redshift bins. In each 
redshift bin we calculate HXLFs both for all spectroscopically identified sources together 
{blue squares; hereafter we will sometimes refer to this sample as either "full" or "all spectral 
types") and for broad-line AGNs alone (r ed diamonds) . Thes e binned HXLFs can be seen 



in Figures [7b— f. They are consistent with iBarger et al.l (l2005l ). Moreover, as in that paper 



the differences between the full and broad-line HXLFs are clear: At high X-ray luminosities 
the full and broad-line HXLFs are consistent, but while the full HXLF continues to rise with 
lower X-ray luminosities, the broad-line HXLF appears to fall off. 
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2 - 8 keV Luminosity {ergs s"') 2 - 8 keV Luminosity (ergss"') 



Fig. 7. — Binned rest-frame 2 — 8 keV luminosity functions for our (a) < 2; < 0.1 local 
sample (see §5.3p and for our distant sample in the low-redsliift bins (b) 0.1 < 2; < 0.4, 
(c) OA < z < 0.8, and (d) 0.8 < 2; < 1.2 and in the high-redshift bins (e) 1.5 < ^ < 3.0 
and (f) 3.0 < z < 5.0. Blue squares (red diamonds) denote the HXLFs for all spectroscopi- 
cally identified sources (broad-line AGNs) in each redshift bin. Error bars indicate the 1 a 
Poissonian uncertainties based on the number of objects in each bin, while arrows denote 
90% (2.3 object) upper hmits. The black dot-dashed curves in each panel (except a) give 
the maximal HXLFs found by assigning the central redshift of that redshift bin to all the 
spectroscopically observed but unidentified sources. The blue solid and red dashed (b— d) 
curves, respectively, show the maximum likelihood fits over 10^^ < Lx < 10^^ ergs s~^ and 
< z < 1.2 using the ILDE model (plotted at the geometric mean of each redshift bin) 
for the full and broad-line data of the distant sample (see §6.2p . The z = extrapolation of 
the fits is shown in (a), and the extrapolation of the full sample fit is also included in each 

ft 
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We also calculate maximal luminosity functions in each redshift bin to estimate the 
effect of incompleteness on our results. For a particular bin, this is achieved by assigning 
all the spectroscopically observed but unidentified sources with valid flux measurements a 
redshift at the center of the bin. This calculation is then performed for each redshift bin, 
giving results which are therefore incompatible but which provide an estimate of the largest 
possible error due to incompleteness. These results can be seen as the black dot-dashed 
curves in Figures [7)d— f. 

Incompleteness will have some effect in each redshift bin but is most likely to bias the 
results for higher redshifts due to the difficulty in getting spectroscopic identifications of low 
luminosity, high-redshift sources. Large incompleteness at thes e redshifts would he lp explain 



the flattening of the HXLFs relative to the low-redshift bins. iBarger et al.l (120051 ) assigned 



photometric re dshifts (where pos sible) to the unidentified sources in the CDF-N and CDF-S 



(working from lZheng et al.ll2004l for the latter) and computed a "spectroscopic plus photo- 
metric" luminosity function. They found the largest discrepancy with their spectroscopic- 
only results to be in the 1.5 < 2 < 3.0 bin, which suggests that most of the unidentified 
sources in those samples lie in this redshift bin. Indeed, a look at Figure H] admits the 
possibility of an underdensity of low luminosity sources in this redshift region due to the 
incompleteness of these samples. 

Another way to understand how our incompleteness could affect our results is to impose 
stricter flux limits such that the remaining sample is ~ 90% spectroscopically complete. We 
performed such a cut to each survey separately and then recalculated our full and broad-line 
AGN HXLFs. There was little variation in either, within the uncertainties, at low redshifts. 
At higher redshifts, the bright end showed no significant change, but the self-imposed limits 
resulted in a complete lack of faint end information. This affirms our conjecture that our 
incompleteness is related to faint, high-redshift sources, but it does not provide us with an 
estimate of the magnitude of the effect. 

While incompleteness is a major source of uncertainty when considering the full HXLF 
at high redshifts, we do not believe the broad-line AGN data suffer from such a problem. 
This is because broad-line AGNs are easier to detect and measure spectroscopically. We 
therefore believe our broad-line AGN data to be fairly complete at least out to z = 3, where 
we have enough data to calculate a reasonably constrained luminosity function. 

An important characteristic of our HXLF for the broad-line AGNs is the significant 
turnover in number density below the luminosity break. A similar turn over in number 



density is not seen the soft X-ray luminosity function of "type-1 AGNs" of iHasinger et al. 



20051 



]). In con t rast to our use of pure spectroscopic redshifts and optical classifications. 



Hasinger et al.l (120051 ) use both spectroscopic and photometric redshifts, and they adopt a 
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mixed classification scheme for "type-1 AGNs" tliat includes any object optically classified 
as a broad-line AGN as well as any object with Lx > 10^^ ergs s~^ and a Chandra - specific 
hardness ratio (Chard — CsnfJ /(Ch ard + Csoft) < —0.2, where Chard and Csoft are the 2 — 8 keV 



and 0.5 — 2 keV count rates (see lSzokoly et al.l 12004 ). (They also employ a flux-dependent 



completeness correction in the calculation of their HXLF, which is known to be biased 
when the remaining sources have not been identified for non-random reasons, though these 
corrections are small in their case.) They adop t this mixed clas sification scheme because 
according to the simple unified AGN model (e.g., Antonuccilll993l ) there should be a one-to- 
one correspondence between the optically unobscured (broad-line) and the X-ray unobscured 
(type-1) AGNs, with a similar relation for the obscured sources. However, studies have 
shown that there i s a mismatch between optical and X-ray identifications of 10 — 20% (e.g., 
Garcet et al.l 120071 and references therein; Trouille et al. 2009, in preparation). The physical 
mechanism for this mismatch is not yet fully understood, and thus we avoid the use of mixed 
classifications, which may introduce unknown biases and effects. 

It has been suggested that pure optical classification schemes may misidentify broad-line 
AGNs as optically "no rmal" galaxies (i .e., our star former and absorber classes), particularly 
at lower luminosities (IGilli et al.l 120071 ). For example, it has been shown that local Seyfert 
2 galaxies may not be prop erly identified as such due to the host gala x y light overwhelming 
that of the nuclear region (IMoran et al.l 120021 : ICardamone et al.l l2007l ) . iBarger et al.l (120051 ) 
tested this galaxy dilution hypothesis using their CDF-N data by comparing the nuclear 
UV magnitudes of the sources with t he 0.5 — 2 keV fluxes, which are known to be strongly 
correlated for broad-line AGNs (e.g.. IZamorani et al.lll98ll ). They found that the optically 



identifled broad-line AGNs showed this correlation, while the other classes did not (their 
nuclei were much weaker relative to their X-ray light than wo uld be expected if they were 
similar to the broa d- line AGNs). Rece ntly, ICowie et al.l (120091 ) combined Galaxy Evolution 
Explorer (GALEX: [Martin et al.l[2005[ ) observations with the CLASXS, CDF-N, and CDF-S 
X-ray samples to determine the ionizing flux from 2; ~ 1 AGNs, and they found that only 
the broad-line AGNs are ionizers; all of the non-broad-line AGNs are UV faint. From these 
two lines of evidence we conclude that we are not misidentifying sources as non-broad-line 
AGNs when they are really broad-line AGNs, as one might have expected to happen if the 
broad lines were not visible spectroscopically due to dilution by the host galaxy. 

We also considered the effect that choosing a different -Zgpiit, the redshift separating our 
use of the hard band and soft band samples, would have on our results. Instead of the 
-2^spiit = 3.0 we have been using, we tried a value of Zsput = 1-5. The consequence of this 
is that both the 1.5 < z < 3.0 and 3 < 2; < 5 redshift bins of our HXLF are determined 
by our soft band sample, rather than just the latter. We flnd that this does not have a 
significant effect on the results of this or other sections, and thus we claim our analysis to 



be independent of the choice of -Sgput- 
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Fig. 8 . — Local binned 14 — 195 keV luminosity function {black squares) and iTueller et al. 



(120081 ) 's analytic approximation {black curve) to this quantity for the |6| > 15° BAT sample. 
The red diamonds show the binned 14 — 195 keV luminosity function for the broad-line 
AGNs, and the blue triangles show the same for the X-ray soft sources (defined as having a 
2 - 8 keV flux to 14 - 195 keV of 0.14). 
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5.3. Results for the Local Sample 



Using the |6| > 15° BAT sample, we constructed binned 14 — 195 keV luminosity 
functions for the full sample, for the broad-line AGNs, and for the X-ray soft sources (which 
we take to be sources for which the ratio of the 2—8 keV flux to the 14—195 keV flux is greater 
than 0.14). In Figure [H] we compare the luminosity functions for the broad-line AGNs [red 
diamonds) and for the X-ray soft sources {blue triangles) wit h the total luminosit y function 
{black squares), as well as with the analytic approximation of iTueller et al.l (120081 ) (see their 
E g. 1 and Table 2; bla ck curve). Our result for the total luminosity function reproduces that 
of 



Tueller et al.l fl2008f ). 



We see that the two selections, optical spectral class and X-ray softness, give very similar 
results. Nearly all of the most luminous sources are X-ray soft broad-line AGNs, and it is 
only below the break in the luminosity function at a luminosity of 10^"^ ergs s~^ (14—195 keV) 
that there are signiflcant numbers of X-ray hard sources. Even here the bulk of the sources 
are X-ray soft broad-l ine AGNs, which su ggests that there are not huge numbers of highly 
absorbed sources. As iTueller et al.l (120081 ) pointed out, the total luminosity function has a 
steep er faint-end s l ope t han the extrapolation of the 2 — 8 keV luminosity function to z = 
from lBarger et al.l (120051 ): howe ver, the faint- e nd slo pe of the X-ray soft sample is flatter and 
hence in better agreement with lBarger et al.l (120051 ) . This is reassuring, since it is the X-ray 
soft sources which will appear in the 2 — 8 keV samples. 

We can directly quantify this by constructing the 2 — 8 keV luminosity function from the 
BAT selected local 2 — 8 keV sample (see §2.21 for a description; this is the sample we have 
been referring to as our local sample). We show this luminosity function in Figure [7^ for all 
Lx > 10^^ ergs s~^ spectroscopically identifled sources {blue squares) and broad-line AGNs 
{red diamonds) in the redshift bin < z < 0.1. The comparison of the binned data with 
the extrapolations of our maximum likelihood flts will be discussed later. In Figure [9] we 
compare the luminosit y function for the full sa r nple ( blue squares) with the local 3 — 20 keV 
luminosity function of ISazonov fc RevnivtsevI (120041 ) from RXTE data (with and without 
their incompleteness correction of 1.4; red dotted curves) after conversion to 2 — 8 keV 
assuming a general photon index F = 1.8. We flnd good agreement. 



In the construction of their HXLF, lUeda et al.l (|2003l ) include a sample of 49 AGNs 
from two HE A 1 surveys which cover a region of Lx — z space very similar to that of our 
local BAT sample. Our chosen flux limit of 1.7 x 10~^^ ergs cm~^ s~^ in the 2 — 8 keV band 
is comparable to that of the HE AO 1 MC-LASS survey and only slightly deeper than the 
~ 2.2 X 10^^^ ergs cm^^ s^^ limit (when converted to 2 — 8 keV) of the HEAO 1 A- 2 survey. 



from w hich 28 of their 49 sources are drawn. These samples are included in the lUeda et al. 
( I2OO3I ) total evolutionary HXLF model flts, which are flt out to 2; = 3 using a total sample 
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of 247 AGNs. We show the z = extrapolation of their LDDE model fit in Figured [green 
dashed curve). The curve agrees with our local HXLF only at the faintest lum inosities and 
otherwise lies above our val ues and the c urves of ISazonov fc RevnivtsevI (120041 ) . This may in 
part be due to the fact that lUeda et al.l (120031 ) fit for an "intrinsic" HXLF before absorption 
effects, but it is also likely related to their limited sample size and to problems associated 
with getting an accurate representation of the local luminosity function while simultaneously 
fitting over an extended redshift region. 
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Fig. 9. — Binned rest-frame 2 — 8 keV luminosity function for our < z < 0.1 local sample 
fs ee ^5.31) (blue squares). The re d dotted curves show the local 3 — 20 keV luminosity function 
of ISazonov &: Revnivtsevi (12004J ) using RXTE data (with and without their incompleteness 
correction of 1.4) after converting to 2 — 8 keV assuming a general photon ind ex F = 1.8, 
while the green dashed curve is the z = Q extrapolation of the LDDE model fit of lUeda et al. 
( 120031 ) . converted from 2 — 10 keV to 2 — 8 keV luminosity assuming F = 1.8. 
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6. Hard X-ray Luminosity Functions: Maximum Likelihood Fits 



It is often of particular use to find an analytic form for a HXLF independent of the 
binning procedure described above. A standard approach is to assume some model with 
unknown parameter dependen ce and use the maxim um likelihood method to estimate the 
values which best fit the data ([Marshall et al.lll983l ). To do this, we define the function 



N(\ogLx,z) = (praodclOog Lx, Z) n(\og Lx, Z) -^(^) 



(12) 



where 0modei is our parameter dependent model for the HXLF. N{logLx, z) is simply the 
number density of AGNs per logarithmic luminosity per redshift predicted by the model. The 
best-fit parameter values in a region of the HXLF defined by logLmin < logi^x < logi^max 
and Zmin < z < z^s.x are those which minimize the function 



-2$: In 



A^(log Lx,i,Zi 



J J N{logLx,z) dlogLxdz_ 



(13) 



whos e sum is over all ob jects in the sample which occupy this region of (logLx,-2) space 
(e.g., iMiyaji et al.l |2000| ). The 1 a errors for the parameters are then estimated as the 
increase/decrease needed to achieve AC = 1 from the minimized value, Cmm, while simulta- 
neously minimizing the function over all other parameters. 

The function C will be independent of any overall normalization of (pmodei- To estimate 
such a parameter, we require that the number of objects predicted by the model in the 
given region be equal to the observed number. The 1 a error will then be given by Poisson 
statistics, with a value ~ A/\/N for normalization A and total number of objects A^. 

Unlike a x^-fitting, the maximum likelihood method does not provide any estimate of 
goodness-of-fit. Therefore, we use Kolrnogoro v-Smirnov (K-S) tests as independent estimates 
of the quality of the fits (jPress et al.l [1992) • Beginning with the distribution defined by 
A^(logLx, ;2), 1-D K-S tests are performed on both the logLx and z distributions (found 
by integrating out one of the components and normalizing) and 2-D tests are performed on 
the full, normalized distribution. It should be noted that these K-S tests are designed to 
give the probability that a sample distribution may have been drawn from a given model 
distribution. This interpretation is not strictly true when any of the model parameters have 
been estimated from the sample distribution, as is the case here, and therefore these tests 
should only be thought of as a rough estimate of the goodness-of-fit. 
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6.1. Models 



For each of the model s we consider, the z = behavior is represented as a double power 



law flPiccinotti et al.lll982l ) 



0) 



dlogL 



A 



X 



71 



+ 



72 



(14) 



The manner in which this form evolves with redshift characterizes each class of models. 
The two most basic involve an evolution of either or A. These models are referred to as 
pure luminosity evolution (PLE) and pure density evolution ( PDE) models, resp e ctively, and 



have been studied in the X-ray samples by va rious authors (IMiyaji et al.l l2000l : lUeda et al. 



2003: iHasinger et aL 
Barger et all J2OO5I I 



2OO5I : iBarger et al.l l2005l ) . We begin by considering the form used by 



which allows for independent evolution of each: 



where 



d$(Lx,z) d$(Lx/eL(^),0) 



dlogL 



X 



dlogL 



X eDiz) 



(15) 



X 



e^(z) = (l + ^)^^ 



'l + z 



,pD 



(16) 
(17) 



This model, which for the purpose of this paper we will refer to as the independent luminosity 
and density evolution (ILDE) model, has six free parameters, four of which determine the 
z = behavior (L*, A, 71,72) and two of which contribute to the redshift evolution of the 
luminosity knee and overall density normalization {pL,pD). 

We also consider a luminosity dependent density evolution (LDDE) model with the 
general evolutionary behavior given by 



d<l>(Lx,z) d$(Lx,0) 



dlogL 



X 



dlogL 



X e{z,Lx) 



(18) 



X 



We actually utilize a slight redefinition of the parameters used in barger et all |2005h . who fit for the 
values of the luminosity "knee", Lq, and the density normalization, oq, at z = 1. We define equivalent z = 
parameters, and A, such that = Lo{1/2)p^ and A = ao(l/2)P^. This has no effect on the actual 
behavior of the model. 
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A useful form of e{z, Lx) introduced by lUeda et al.l (120031 ) is 



e(z, L 



X) 



(1 + z < z,{Lx) 

e{z,)[{l + z)/{l + z,)]P^ z>z,{Lx) 



(19) 



z* 



Lx<L, 



(20) 



This model has nine parameters: the four z = parameters and five others governing the 
redshift evolution {pl,p2, z*, La,a). It is usually the case that pi > and p2 < 0, allowing 
for a density increasing as (1 + z^^ up to so me luminosi t y dep endent cut-o ff, Zr.(Lx), after 
which the density drops as (1 + 2;)^^. Unlike lUeda et al.l (120031 ) and others (IHasinger et al. 



ps 



2OO5I : [Silverman et al.l 120081 ). we do not fix any of these evolution parameters at a constant 
value. While this complicates the procedure slightly, it allows more fiexibility in fitting the 
data. 



6.2. Results for the Distant Sample 

In Table [3] we give the maximum likelihood fit parameters obtained using the ILDE 
model (which we only fit over < z < 1.2) for all spectroscopically identified sources and 
for broad- line AGNs. We do the fits over the luminosity int erval 10"^^ < Lv < 1 0^^ ergs s^^. 



The parameters for all spectral types are consistent with iBarger et al.l (120051 ). given the 
uncertainties. While this is true of both the z = and evolution parameters, it is worth 
noting that the values of pD = — 0.94lQ;y5 for the full sample and pD = — 0.8ltQ;g2 for the 
broad-lines do suggest some density evolution. We plot the fits (all spectral types - blue solid; 
broad-line AGNs - red dashed) in Figures [7)d— d and the z = extrapolations in Figure [7^ 
(note that these fits do not include the local data shown in (a)). We also include the z = 
extrapolation of the full sample fit in each of Figures [7)d— d (black dotted). Visually, the 
curves provide a plausible, albeit simplified, match to the binned data, including the local 
luminosity function. While the K-S tests indicate that the ILDE model may not adequately 
fit the unbinned data for all spectral types, the fit to the unbinned data for broad-line AGNs 
is much better, with 1-D K-S tests on the luminosity and redshift distributions at 91% and 
81%, respectively, and a 2-D test at 57%. 

In Table H] we give the maximum likelihood fit parameters obtained using the LDDE 
model for all spectroscopically identified sources and for broad-line AGNs. In each of the two 
cases we have fitted over three redshift intervals (0 < 2; < 1.2, < 2; < 3, and < z < 5). 
We give the results for each interval in the table. While the parameters of the < 2; < 3 and 
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< 2; < 5 fits to all spectral types are generally consistent with each other, the parameters 
of the < z < 1.2 fit to all spectral types are typically quite different, reflecting the strong 
negative density evolution that occurs at higher redshifts. In fact, it is the value of p2, which 
characterizes the high-redshift density evolution, that most distinguishes each of the three 
flts, showing a continual decrease as higher redshift data are included. 

Of our two high-redshift flts, the < 2; < 3 flt to all sp ectral types shows the most consis 



tency with other HXLF studies that use an LDDE model jUeda et al.ll2003l : iLa Franca et al. 



2005: Ebrero et al 



2008 



Silverman et al 



studies, as lUeda et al. J2003h. lEbrero et al.l (120081 ) , and [Silverman et al.l (120081 ) have flts out 
to z = 3, while iLa Franca et al.l (120051 ) have flts out to ^ = 4. We show the most consistency 
with [Silverman et al.l (120081 ) . signiflcantly differing only in the overall normalization A. Our 
value of logy4 = —6.18 ± 0.04 is noticeably lower than their log A = —6.08 ± 0.02. 



This re dshift interval is c ompa rable to those 



It is important to note that, like [Silverman et al.l (120081 ). we do not correct for ab - 
sorption in our HXLF mode l flts. This i s in co ntrast to the works of lUeda et al.l (120031 ). 
La Franca et al.l (120051 ) , and lEbrero et al.l (120081 ) , all of whom flt for the effects of absorp- 
tion separately and determine "intrinsic" , unabsorbed luminosity functions. In spite of this, 
we typically show good agreement with most of their best-flt parameters, given the uncer- 
tainties. Howev e r, the re are some notable exceptions. Similar to our disagreement with 



Silverman et al.l (120081 ) . our estimate of A is lower than each of their values. We flnd our 
value of log L^, = 44.28lo;28 (in units of ergs s~^) to be larger than that0 found by Ebrero et al. 
(I2OO8I), logL, = 43.8 3 +°-°^, by ~ 2a, though consist ent with log L, = 44.17 ± 0.18 of Model 



43.86^°;^^, 



4 of iLa Franca et al.l (120051 ) and nearly with that of lUeda et al.l (120031 ) . logL* 
given the uncertainties on these values. Also, our v alue of logLp = 4 4.68j^ n'i9 is nearly 
consis tent with the flxed value l o gL„ = 44.52 used by lUeda et al.l (120031 ) and lEbrero et al. 



(120081 ) , though iLa Franca et al.l (120051 ) flt this parameter freely and flnd a larger value of 
logL, = 45.681° js. 

Of interest is how any disagreements in best-flt parameters are manifest in the z = 
1.5 — 3.0 region of the HXLF, particularly at the faint end where it is difficult to make 
accurate redshift measurements. In Figure [10] we show our binned HXLF in this redshift 
bin {blue squares), along with our best-fit < z < 3 LDDE model {blue solid curve). 
This is pl otted at the geometric mean of the redshift bi n, along with the best- fit LDDE 



models of Ueda et al. 



dashed), [Ebrero et al 



(I2OO3I : green short dashed- dotted), [La Franca et al.l (l2005l. r ed long 
(I2OO8I : black long dashed- dotted), and [Silverman et al.l (l2008l : purple 



"^W e hav e converted a l l of th e 2 — 10 keV luminosity parameters of lUeda et al.l (|2003l ). [La Franca et al 
|2005l ). and [Ebrero et all |2008h to 2 - 8 keV luminosities assuming T = 1.8. 
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sh ort dashed). As descri bed above, the shape of our HXLF agrees very well with that 
of ISilverman et al.l (I20081) but dif f ers in the overall normalization. Our fit also matches 
very well with that of lUeda et al.l (120031 ) at the bright end of the luminosity function, but 
their fit lies above ours below the lu minosity break. This d iscrep ancy at the faint en d is 
even more apparent in the curves of iLa Franca et al.l (120051 ) and lEbrero et al.l (120081 ) . It 
is also clear from Figure [10] that the LDDE fits of each of these studies lies above our 
binned determinations of the HXLF at the faint end, while showing some agreement at 
the bright end. The incompleteness of our sample at faint fluxes is likely a large fa c tor in 
these differences, as discussed in §5.21 Our closer agreement with [Silverman et al.l (120081 ) 
suggests that some discrepancy with the other studies may also come from trying to directly 
compare our HXLFs with their absorption corrected, intrinsic HXLFs. This can also be 
seen by comparing the curves to the maximal binned HXLF in the z = 1.5 — 3.0 bin, 
found by assigning all the spectroscopically observed but unidentified sources with valid 
flux measurements a redshift at the ce nter of the bin. This is gi ven by the black do tted 
curve in Figure [TOl While the HXLFs of iLa Franca et al.l (I20051 ) and lEbrero et al.l (I20081 ) are 
similar to the maximal HXLF at l ow luminosities, t hey d o no t match the overal l shape for 
all luminosities, and the HXLFs of [Silverman et al.l (l2008l ) and lUeda et al.l (120031 ) appear to 
agree with the maximal curve only at higher luminosities. 
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Fig. 10. — The squares give the binned rest-frame 2 — 8 keV luminosity function in the 
redshift region 1.5 < z < 3.0. Error bars indicate la Poissonian uncertainties on the number 
of objects in each bin. The blue sohd hne is the best-fit LDDE model of our distant sample 
over the redshift region < 2; < 3, plotted at the geometric mean of the redshift bin, and 
the black dotted line gives the maximal binned HXLF for this redshift region (found by 
assigning the central redshift of the bin to all unidentified sources in the sample). The 
other curves give the best-fit LDDE models of other studies, corrected from 2 — 10 keV to 

1.8 and plott ed at the geometric ra ean of the redshift 

red long dashed; 



bin ( 


Ueda et al. 


2003 


Ebrero et al. 


200 


8 - b 



- black long dashed- dotted: ISilverman et al.ll2008l - purple short dashed). 
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In addition to their fit over 0.2 < z < 3.0, ISilverman et al.l (120081 ) find best-fit LDDE 



parameters over 0.2 < z < 5.5. Their sample contains 31 AGNs wi th z > 2>, while our soft 



band sample has 41 AGNs with 2; > 3 (2 of which lie above z = b). ISilverman et al.l (120081 ) 
choose to fix all but 4 of their parameters (with the others set to the values found by their 
0.2 < 2; < 3.0 fit), while all of the parameters in our < 2; < 5 fits are allowed to vary 
freely. This makes a direct comparison somewhat difficult. However, they find a value of A 
which is in very good agreement with our own. They also find that a steeper value of p2 is 
necessary to match the high-redshift data, with their estimate of p2 = — 3.27;1]q 34 showing 
consistency with our p2 = — 2.83lo24 within the uncertainties. The only major disagreement 
between our fits is in the value of a. Their a = 0.333 ± 0.013 is much larger than our 
a = 0.208lo'o}9. This difference is likely caused by their fixing of La, a very complementary 
parameter. Despite this, comparisons of the two LDDE models in the redshift bin 3 < z < 5 
show good agreement. 

To give a sense of the full redshift evolution of our fits, we plot our < 2; < 5 LDDE 
model fit in different redshift bins in Figure [11] (all spectral types - blue solid; broad-line 
AGNs - red dashed). We show the z = extrapolations in Figure [TTk. We also include the 
z = extrapolation of the full sample fit in each of Figures [TTb— f {black dotted). (We note 
that our LDDE model fits over other redshift intervals show similar features.) The fits to all 
spectral types show consistency with the binned luminosity functions over the full redshift 
interval fitted, but the z = extrapolations are in relatively poor agreement with the local 
luminosity function in Figure [TTh . The LDDE model fits in Figure [TT] reaffirm the notion 
that the broad-line AGN population density drops at low luminosities relative to that of the 
full population, particularly at higher redshifts. 
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Fig. 11. — Binned rest-frame 2 — 8 keV luminosity functions for our (a) < 2; < 0.1 local 
sample (see §5.3p and for our distant sample in the low-redsliift bins (b) 0.1 < z < 0.4, 
(c) OA < z < 0.8, and (d) 0.8 < 2; < 1.2 and in the high-redshift bins (e) 1.5 < ^ < 3.0 
and (f) 3.0 < z < 5.0. Blue squares (red diamonds) denote the HXLFs for all spectroscopi- 
cally identified sources (broad-line AGNs) in each redshift bin. Error bars indicate the 1 a 
Poissonian uncertainties based on the number of objects in each bin, while arrows denote 
90% (2.3 object) upper hmits. The black dot-dashed curves in each panel (except a) give 
the maximal HXLFs found by assigning the central redshift of that redshift bin to all the 
spectroscopically observed but unidentified sources. The blue solid and red dashed curves, 
respectively, show the maximum likelihood fits over 10^^ < Lx < 10"^^ ergs s~^ and < z < 5 
using the LDDE model (plotted at the geometric mean of each redshift bin) for the full 
and broad-line data of the distant sample. The z = extrapolation of the fits is shown in 
(a), and the extrapolation of the full sample fit is also included in each higher redshift panel 
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Fig. 12. — Binned rest-frame 2 — 8 keV luminosity functions for our (a) < 2; < 0.1 local 
sample (see §5.3p and for our distant sample in the low-redshift bins (b) 0.1 < z < 0.4, (c) 
0.4 < z < 0.8, and (d) 0.8 < z <\.2. Blue squares denote the HXLFs for all spectroscopically 
identified sources in each redshift bin. Error bars indicate the 1 a Poissonian uncertainties 
based on the number of objects in each bin, while arrows denote 90% (2.3 object) upper 
limits. The red solid curve gives the fit to the distant sample over < z < 1.2 using the 
ILDE model, while the green long-dashed curve uses the LDDE model (the curves are shown 
at the geometric mean of each redshift bin). The z = extrapolation of the fits is shown in 
(a). Also in (a), the red dotted curves at the high luminosity end and the green short-dashed 
curves at the low luminosity end show the 1 a spreads for the z = extrapolations using 
the ILDE model and the LDDE model, respectively. 
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While the fits to all spectral types are likely biased at high redshifts due to potentially 
large incompleteness, the broad-line AGNs are relatively free from such effects. The LDDE 
model fits to the broad-line data over < z < 3 and < 2; < 5 show good agreement 
with the binned estimates, as seen in Figure [TT] for the < z < 5 case, and the K-S tests 
suggest very good fits to the unbinned data. The K-S tests are also favorable for the fit to 
the broad-line data over < z < 1.2. However, the best-fit values of z*, La, and a for the 
< 2; < 1.2 fit are quite different than those for the two high-redshift fits, and this results in 
curves which are rather cuspy and unlikely to be a good description of the actual population 
of low-redshift broad-line AGNs. In general, we take the LDDE model fit over < z < 3 to 
be our best parameter set for broad-line AGNs in the distant sample. 

Figures [T2b—d compare the fits over < 2; < 1.2 of the ILDE {red, solid curve) and the 
LDDE {green, long-dashed curve) models to all spectral types for a finer binned HXLF. The 
models agree very well in the z = 0.8 — 1.2 bin but begin to diverge in the lower redshift 
bins. This can also be seen in Figure [T2h . where we compare the z = extrapolations 
of the fits (and their 1 a spread) with our local HXLF. Here it can also be seen that the 
extrapolation of the fit using the ILDE model matches the HXLF better at low luminosities, 
while the extrapolation of the fit using the LDDE model matches the HXLF better at high 
luminosities. 



6.3. Results for the Distant Plus Local Sample 

We repeated our analysis on the combined distant plus local sample. In Table [5] we 
give the maximum likelihood fit parameters obtained using the ILDE model (which we again 
only fit over < 2; < 1.2) for all spectroscopically identified sources and for broad-line 
AGNs. We do the fits over the luminosity interval 10^^ < Lx < 10^^ ergs s~^. The fit to all 
spectral types shows virtually no deviation from the fit in §6.21 other than a reduction in the 
uncertainties of some of the parameters. There is some variation in the parameters for the 
fit to the broad-line AGNs compared to the fit in §6.21 but, with the exception of a slight 
increase in 72, these are not statistically significant. 

In Table [6] we give the maximum likelihood fit parameters obtained using the LDDE 
model for all spectroscopically identified sources and for broad-line AGNs. In each of the two 
cases we have fitted over three redshift intervals {0 < z < 1.2, < 2; < 3, and < 2; < 5). We 
give the results for each interval in Table [61 For the fits of all spectral types, the parameters 
in Table [6] are generally consistent with the parameters in Table H] but suggest a lower value 
of logL* and larger log A and pi. The changes in the former two parameters are the result 
of bringing the HXLF into better consistency with the local data, while the latter corrects 
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for the change this would bring about at higher redshifts. 

There is also a general trend among the Table M fits to broad-line AGNs for larger 
values of the faint and bright end slopes (governed by 71 and 72). This would be a mild 
improvement over the Table H] fits as it would bring the bright ends of the fits to the broad- 
line AGNs into greater consistency with the corresponding fits to the full sample. The faint 
ends of the fits to the broad-line AGNs, though potentially steeper, still lie well below the 
full population. The < 2; < 1.2 fit, already problematic when fitting the distant sample, is 
further complicated by the addition of the local data. The z = parameters show a strong 
break with those of the high-redshift fits and the K-S tests indicate an unsatisfactory fit. 
Fortunately, the < ;z < 3 and < z < 5 fits provide a better description, with (logLx, z, 
2-D) K-S probabilities in the < ^ < 1.2 region of (0.97,0.59,0.54) and (0.99,0.56,0.49), 
respectively. 

In Figure [13] we compare the fits over < 2; < 1.2 for the distant plus local sample using 
the ILDE {red solid curve) and LDDE {green long-dashed curve) models with the binned 
local HXLF. Each model now agrees well with the local data at the faint end. Additionally, 
we show the 1 cr spreads on the bright end slope to demonstrate the consistency with the 
binned data in that region as well. At higher redshifts, the LDDE model fits are very similar 
to those shown in Figure [TTl 

Combining the local sample with the distant sample generally improves the quality of 
the LDDE model fits over all redshift ranges. We therefore recommend using the LDDE fits 
of the combined sample over < 2; < 3 for all the spectroscopically identified sources and 
for the broad-line AGNs as our best overall fits. The LDDE fits over < z < 5 are also 
quite reasonable and could be useful when a description to very large redshifts is needed. 
However, it must be noted that when dealing with the sample of all identified sources, the 
fits may be biased by incompleteness, with the greatest effects occurring for z > 1.5. For 
this reason, if one is only interested in redshifts z < 1.2, then the LDDE fit of the combined 
sample over < z < 1.2 for all identified sources (but not for broad-line AGNs) is likely to 
provide a better description of the data. 
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Fig. 13. — Binned rest-frame 2 — 8 keV luminosity function for our < z < 0.1 local sample 
(see §5.31) {blue squares). Error bars indicate the 1 cr Poissonian uncertainties based on the 
number of objects in each bin, while arrows denote 90% (2.3 object) upper limits. The 
red solid curve gives the z = extrapolation of the fit to the distant plus local sample over 
< 2; < 1.2 using the ILDE model, while the green long-dashed curve uses the LDDE model. 
The red dotted curves and the green short-dashed curves show at the high luminosity end 
the 1 a spreads for the z = extrapolations using the ILDE model and the LDDE model, 
respectively. 
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6.4. Space Density Evolution 



In addition to tlie luminosity function, the comoving space density of AGNs as a function 
of redshift can be useful in understanding the evolution of both the full and broad-line AGN 
populations. These can be derived from an unbinned HXLF (which is simply the comoving 
space density per unit logarithmic luminosity) by integrating over the luminosity region of 
interest at a particular redshift. In principle, a binned estimate can also be found by first 
computing a binned HXLF and multiplying by the size of the luminosity bin. We will not 
follow this approach here, as our method of estimating the binned HXLF as described in 
^ assumes that it varies little over the bin, and we will be choosing luminosity bins large 
enough such that this assumption breaks down. 

We will th erefore estimate the comov i ng space density acc ording to the traditional 
1/Kiax method JSchmidtl EoGsl : lFeltenlll976l : lAvni fc Bahcalllll98o[ ). At a redshift z in a bin 
Zi < z < Z2 and in the luminosity range log Li < log Lx < log L2 this method gives 



(21) 



where the sum is over all objects in this region of {logLx, z) space and 



Knax,i = j ni\0gL^,,,z') — {z')dz' 



(22) 



for a source with luminosity logLx.i. The factor Knax,i is the maximum accessible volume of 
each source as dictated by its luminosity and the survey details incorporated in f2(logLx, z) 
(see §l]for details). The estimate of p{z) is then just the sum of the individual contributions 
to the total space density. 

In Figures [T^ and we plot our binned estimates for the comoving space density 
for the full and broad-line populations, respectively, in four luminosity regions. The er- 
rors are estimated as the la Poissonian uncertainties based on the number of objects in 
each bin. Also plotted are the curves derived from integrating the unbinned HXLF over 
the corresponding luminosity intervals. We use the LDDE model with parameters deter- 
mined by the fit of the local plus distant sample over 2; = — 5, as given in Table [61 
For both populations it is clear that the peak density of higher luminosity sources occurs 



suDDorts the idea of "cosmic downsizing" (ICowie et al. 
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broad-line best-fit curves displayed separately in Figures [T^ and UM i is given in Figure [151 
It is clear that while the broad-line population makes up a significant portion of the full 
population at higher luminosities for all redshifts, it is only below the redshifts at which the 
density peaks that faint broad-line sources make a sizable contribution. This is similar to 
what is seen in the redshift evolution of the luminosity functions. 

A close look at the binned data for the lowest luminosity region in Figure UM ) also 
reveals an interesting trend. Unlike every other luminosity class for both the full and broad- 
line AGN samples, this class of broad-line AGNs is only now reaching its peak densitjl^l. This 
is interesting to note in the context of the binned HXLF in Figure [111 in which it appears 
the low luminosity broad-line AGNs in the local sample have a noticeably larger density 
than in any of the higher redshift bins. It is clear from Figure [Hb, though, that the local 
data simply continue the upward trend in density and are not, in fact, atypical of the data 
as a whole. As discussed in §5.21 we believe our broad-line identifications are robust and not 
significantly affected by host galaxy dilution, and that this is therefore an actual physical 
trend in the low luminosity broad-line AGN data. 



^It should be noted that there is some disagreement between the binned data and the model fit in this 
luminosity region. The discrepancy disappears, however, if the best-fit parameter values for the z = — 1.2 
fit over the distant plus local sample are used, rather than the z = — 5 fit shown, as the z = — 1.2 fit is 
the most sensitive to the low-redshift data. 



-43- 




O 42<logLx<43 

O 43<logLx<44 

□ 44 < log Lx < 45 

A 45 < log Lx < 46 



J I I L_ 



_I I I L_ 



2 3 
Redshift 




Fig. 14. — Comoving space density in the given luminosity intervals for (a) all spectral 
types and (b) broad-line AGN only. The error bars indicate the la Poissonian uncertainties 
based on the number of objects in each bin. The curves give the LDDE model fit of the 
HXLF (using the distant plus local sample over z = — 5 for the appropriate population) 
integrated over the intervals log Lx = 42 — 43 (red dotted) , log Lx = 43 — 44 (blue dashed) , 
logLx = 44 — 45 (green dashed-dotted), and logLx = 45 — 46 (sohd black). 
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Fig. 15. — The upper (lower) curve of each pair gives the co moving space density for the 
full (broad-line) AGN population for each luminosity interval. They are determined by 
integrating the LDDE model fit of the HXLF (using the distant plus local sample over 
z = — 5) over the intervals log Lx = 42 — 43 (red dotted), log Lx = 43 — 44 (blue dashed), 
logLx = 44 — 45 (green dashed-dotted), and logLx = 45 — 46 (sohd black). 
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7. Summary 



In this paper we presented rest-frame hard (2 — 8 keV) X-ray luminosity functions 
constructed from a distant sample drawn from five highly spectroscopically complete surveys: 
CDF-N, CDF-S, CLASXS, CLANS, and ASCA. Three of thes e surveys (CDF-N, CLASXS, 
and CLANS) compose our OPTX project JXrouille et a llboo pl). The addit ion of the CLANS 
field substantially incr eases the sample o yer th e one used by iBarger et al.l (120051 ) . Following 
the general method of iPage fc Carreral (120001 ) . we calculated binned HXLFs out to ^ ~ 5 
for all the spectroscopically identified sources and for the broad-line AGNs. We found that 
incompleteness is still a major source of uncertainty for the full (but not for the broad-line) 
HXLFs at 2; > 1.5, potentially generating flatter HXLFs for lower luminosities than in reality. 

A primary goal of our paper was to fit these data with evolutionary models. However, 
because fits using only distant data may extrapolate poorly to z = 0, we also constructed a 

ected from the very hard (14 - 195 keV) SWIFT BAT all-sky X-ray 



local 2 - 8 keV HXLF se 



survey of iTueller et al.l (120081 ) so we could see how the inclusion of local data changes the 



model fits. We used the maximum likelihood method with the combined distant plus local 
sample, as well as with the distant sample alone, to estimate parameters for two different 
types of models, one which allows for separate luminosity and density evolution (the ILDE 
model) and the other which provides for a luminosity dependent density evolution (LDDE). 

We found that the best-fit ILDE parameters over < 2; < 1.2 for the distant sam- 
ple alone suggest some density evolution, particularly for the broad-line AGN population. 
The best-fit LDDE parameters over < z < 1.2 also give a reasonable fit, further sup- 
porting density evolution at these low redshifts. With the exception of the normalization 
A, our LDDE f i t para meters over < 2; < 3 show very good agreement with those of 
(I2OO8I). Ther e is also fairly good ag reeme nt with the " i ntrins ic" HXLF 



Silverman et al. 



functions of lUeda et al.l (120031 ) , iLa Franca et al.l (120051 ) , and lEbrero et al.l (120081 ) , with the 
exception of some discrepancies in the values of L^., A, and La- This causes some disagree- 
ment at the faint end of the HXLF for z > 1.5, which is likely due to potentially significant 
incompleteness effects caused by the difficulty in spectroscopically observing high-redshift, 
low luminosity sources. Fortunately, as we believe our broad-line AGN sample is relatively 
free of these complications and therefore fairly complete, we have best-fit LDDE parameters 
which provide a very good match to the broad-line AGN data up to 2; ~ 3 and a reasonable 
fit out to 2; ~ 5, as well. 

The fits using the ILDE and LDDE models for the distant plus local sample are typically 
quite similar to the fits without the local data when considering the uncertainty on the 
parameters, but in some cases they do suggest some slight modifications. The fits using the 
LDDE model for the spectroscopically identified sample over < 2; < 3 and < 2; < 5 
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are consistent with slightly lower values of the luminosity break and with larger values 
of the normalization A and low-redshift density evolution parameter pi. This brings the 
z = extrapolations into better agreement with the local binned data while maintaining 
the congruity at high redshifts. Also, in all cases the fits of the broad-line AGN data show 
a mild increase in the values of the faint and bright end slopes, 71 and 72, which brings the 
bright ends of these fits into greater consistency with the fits of the full data while retaining 
the discrepancy at the faint end. We beheve these to be improvements over the fits of the 
distant sample alone and therefore recommend using the LDDE parameters of the combined 
sample over < 2; < 3 for all the spectroscopically identified sources and for the broad-line 
AGNs as our best overall fits, while still noting the possible effects of incompleteness on the 
fits of the full sample. 

When examining both the binned an unbinned HXLFs, constructed with and without 
the local data, we confirm that at all redshifts the population of broad-line AGNs follows a 
fundamentally different form than that of the full population. While the broad-line AGNs 
account for a majority of the overall AGN space density at high luminosities, their density 
drops significantly for lower luminosities. 
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Lucile Packard Foundation (A. J. B.). L. T. was supported by a National Science Foun- 
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Table 3. Maximum Likelihood Fit Parameters for All Spectroscopically Identified Sources 

and for Broad-Line AGNs Using a ILDE Model 



ILDE 




All Identified 


BLAGN 




(0 < ^ < 1.2) 


(0 < ^ < 1.2) 


logL/ 
logA^ 

71 

72 


42 qi+0-25 

^Z.yi_Q 25 

-4.319t°;°^o^ 

0.498tH^^ 
2 in+o-25 

^•^^'-'-0.23 


42.63l°|2 
-4.340tH^i 

-i.osl^i? 


pL 
pD 


4 n4+"-'^5 

^•'-'^-0.86 


3.54t°:^^ 

-o.8ii°:i 


Pks 
(L,, z, 2D) 


(0.93,0.003,0.05) 


(0.91,0.81,0.57) 



^In units of ergs s ^ 
'^In units of Mpc^^ 



Table 4. Maximum Likelihood Fit Parameters For All Spectroscopically Identified Sources and for Broad-Line AGNs 

Using a LDDE Model 



LDDE 




All Identified 


BLAGN 




{0<z < 1.2) 


(0 < z < 3) 


(0 < z < 5) 


(0 < 2 < 1.2) 


(0 < 2 < 3) 


(0 < 2 < 5) 


7i 

72 


AO oq+0.36 

-5 846+°-°^'^ 

0. 540_o.o50 
'-'■^^-0.16 

1 ofi+o.ao 

1. yo_Q 25 


-fi 1 76+"°37 

n 047+0.066 

9 91 +0.21 
'^■'^-'--0.21 


44 4n+0-13 

-6. 1401^1 
0.8721°:°^° 

9 qc+0.24 


44.Uy_i 34 
f- Q99+0.091 

^•^^^-0.101 

n 71 +°-^^ 

•J-' ^-0.62 
1 79+0.98 
-^■'^-0.50 


4q qi+o.20 

O.UiO_o 063 

0.25l°i? 
1.841°:}? 


4'? 81 +0.22 

_c: c:in+0.056 
0.diU_o.060 
n 14+0.190 
'-'■-'-^-0.19 
1 77+0.17 
'-0.15 


pi 
p2 

z* 
logL„^ 

a 


5.84lo;g9 

'J-'^O^Q 52 

n 047+0.055 
43 921 

n QQ7+0.069 
U.OO / _o,o87 


^■^^-0.37 

-1.70l°J^ 

1-801°:}^ 

44.681°;}^ 


q f;i+0.49 
"^■"-'--0.33 
9 QO+0.23 
Z.OO_Q 24 

9 1 0+0.55 
^- J^O-0.29 

45.09+g-f7 
0.208lHi^ 


q qi +0.85 
"^■"^-'--0.87 

_4 0+2.0 

1 099+0.044 
i.UZZ_o 

43.511°-?^ 
92+°-^^ 

*-'-^^-0.22 


3.i6i°:i; 

_q 1 c:+1.41 
•J- -'-^-0.70 

^•^^-0.46 

44.571°-}° 
0.422l°:°«? 


q 90+0.40 
"^■^'-'-0.49 

_q 4c:+0.41 
<5-40_o 41 

2.241°:- 

44.5391°:°^^ 

S90+°-°''° 
u.oyu_o,o33 


Pks 
{L,,z,2D) 


(0.99,0.04,0.23) 


(0.93,0.04,0.08) 


(0.98,0.10,0.18) 


(0.91,0.87,0.95) 


(0.97,0.96,0.96) 


(0.98,0.59,0.43) 



*In units of ergs s ^ 
'^In units of Mpc~^ 
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Table 5. Maximum Likelihood Fit Parameters for All Spectroscopically Identified Sources 
and for Broad-Line AGNs Using a ILDE Model (Distant Plus Local Sample) 



ILDE 




All Identified 


BLAGN 




(0 < ^ < 1.2) 


{0<z< 1.2) 


logL/ 
log 

71 

72 


42 92+°-^^ 

A Q9O+0.046 

0.5001°:°^ 

2.081°;- 


42 70+°-is 
-4 267+°-°^^ 

-0.091 

-0 67+°-^^ 

-0.33 
-'-•yU-0.12 


pL 
pD 


o qs+0.47 
'-'■^-'--0.52 


4 04+0-43 

^••^^-0.22 
-■-■"^^-0.53 


Pks 

(iv,..:.2D) 


(0.77.0.005.0.09) 


(0.13.0.90.0.11) 



^In units of ergs s 
'^In units of Mpc~"^ 



Table 6. Maximum Likelihood Fit Parameters For All Identified Spectral Types and for Broad-Line AGNs Using a 

LDDE Model (Distant Plus Local Sample) 



LDDE 




All Identified 


BLAGN 




(0 < ^ < 1.2) 


(0 < ^ < 3) 


(0 < ;2 < 5) 


(0 < ^ < 1.2) 


(0 < ;2 < 3) 


(0 < ;2 < 5) 


log L^"^ 

log /I 

7i 

72 


4.05Z_Q Q4g 

0.685lo:o48 
1.9361H?^ 


A 'X 00+0-22 
-fi 060+° °^^ 

U.UUU_Q 
■^■-^^-0.20 


44.078_o.o97 

-6.117^Hi 
0.9561°:°^^ 

2.301°:?« 


A A fi/l+0-31 

44.b4_o.2i 

-7 1 89+°-°^''^ 
' --"-o^-o.ogi 

i.iyb_o.o64 

q 1 99+1.40 
O.izz_o 55 


A'X fi/l+0-18 

43.b4_o.i9 

O.O^O_0 060 

0.451}^ 
2.15lS:it 


A'X c;/l+0-24 

43.54_o.i3 

V. /i9c:+0.055 
d.^ZO_o.o58 

0.39l°:J« 
2.091°:}^ 


pi 
p2 

z* 

logLa^ 

a 


6.42t°:i 

1 1 1 +0.39 
-'-■-'--'--0.34 

n qc;Q+o.049 

43.887;^^? 
0.6941°:^^^ 


r c:q+0.41 
^■^°-0.52 
_1 34+0.29 
-'-■"^^-0.24 

1-691^:}^ 

44.68+°;?^ 

n oQo+0.049 


r 1 q+0.35 

^■^"5-0.21 

9 r -1+0.25 
^■'^-'--0.16 
1 064+0-198 

44.889+°;?i5 
255+°-°^^ 

U.Zdd_o 022 


4.801°:!^ 

u.yuu_o.o30 
43.434l°:^26 

1 9O+0.14 
-•-■^^-0.21 


4.93tgi? 
-2 lfi+°-^3 

^•-■-"-0.43 
9 9q+0.11 
^-^"-"-O.l? 

44.566;n47 
n c,c,q+o.65o 

U.ddO_o.o26 


r r)9+0.43 
^■'-'^-0.31 

-2 70+0-28 

i.96;°i^ 

44.4781^^^ 

n c,qc,+0.023 
U.dOd_o.o77 


Pks 


(0.97,0.02,0.14) 


(0.91,0.07,0.15) 


(0.97,0.13,0.43) 


(0.69,0.77,0.31) 


(0.99,0.86,0.85) 


(0.97,0.44,0.40) 



^In units of ergs s 
'^In units of Mpc~^ 
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